A cost‐effective curvature calculation approach for interfacial flows on unstructured meshes

We present a simple and cost‐effective curvature calculation approach for simulations of interfacial flows on structured and unstructured grids. The interface is defined using volume fractions, and the interface curvature is obtained as a function of the gradients of volume fractions. The gradient computation is based on a recently proposed gradient recovery method that mimicks the least squares approach without the need to solve a system of equations and is quite easy to implement on arbitrary polygonal meshes. The resulting interface curvature is used in a continuum surface force formulation within the framework of a well‐balanced finite‐volume algorithm to simulate multiphase flows dominated by surface tension. We show that the proposed curvature calculation is at least as accurate as some of the existing approaches on unstructured meshes while being straightforward to implement on any mesh topology. Numerical investigations also show that spurious currents in stationary problems that are dependent on the curvature calculation methodology are also acceptably low using the proposed approach. Studies on capillary waves and rising bubbles in viscous flows lend credence to the ability of the proposed method as an inexpensive, robust, and reasonably accurate approach for curvature calculation and numerical simulation of multiphase flows.     

Gibbsian dynamics and ergodicity of magnetohydrodynamics and related systems forced by random noise

The magnetohydrodynamics system consists of the Navier-Stokes equations from fluid mechanics, coupled with the Maxwell’s equations from electromagnetism through multiples of non-linear terms involving derivatives. Following the approach of [1 Weinan, E., Mattingly, J. C., Sinai, Y. (2001). Gibbsian dynamics and ergodicity for the stochastically forced Navier-Stokes equation. Commun. Math. Phys. 224:83–106. DOI:10.1007/s002201224083.[Crossref], [Web of Science ®] , [Google Scholar]], we prove the existence of a unique invariant measure in case the forcing terms consist of the cylindrical Wiener processes with only low modes. Its proof requires taking advantage of the structure of the non-linear terms carefully and is extended to various other related models such as the magnetohydrodynamics-Boussinesq system from fluid mechanics in atmosphere and oceans, as well as the magneto-micropolar fluid system from the theory of microfluids.     

Law of the iterated logarithm for random graphs

A milestone in probability theory is the law of the iterated logarithm (LIL), proved by Khinchin and independently by Kolmogorov in the 1920s, which asserts that for iid random variables with mean 0 and variance 1 In this paper we prove that LIL holds for various functionals of random graphs and hypergraphs models. We first prove LIL for the number of copies of a fixed subgraph H. Two harder results concern the number of global objects: perfect matchings and Hamiltonian cycles. The main new ingredient in these results is a large deviation bound, which may be of independent interest. For random k‐uniform hypergraphs, we obtain the Central Limit Theorem and LIL for the number of Hamilton cycles.     

基于有限体积虚拟区域方法的两相流动直接模拟

为提高计算效率，提出有限体积法离散下的虚拟区域颗粒两相流动直接模拟方法.在控制方程中加入相应的虚拟区域源项，保证了颗粒内部的刚体运动特性.该源项中含有颗粒信息部分及流体信息部分.在每次迭代后，对源项中的流体信息部分进行更新，从而更好地保证颗粒内速度的刚体分布.计算静止颗粒圆柱绕流及单个颗粒的沉降过程，验证了算法的准确性.     

Helicity and other conservation laws in perfect fluid motion

In this review paper, we discuss helicity from a geometrical point of view and see how it applies to the motion of a perfect fluid. We discuss its relation with the Hamiltonian structure, and then its extension to arbitrary space dimensions. We also comment about the existence of additional conservation laws for the Euler equation, and its unlikely integrability in Liouville's sense.     

Highly edge-connected regular graphs without large factorizable subgraphs

We construct highly edge-connected $r$-regular graphs of even order which do not contain $r?2$ pairwise disjoint perfect matchings. When $r$ is a multiple of 4, the result solves a problem of Thomassen [4].     

The application of the Lattice Boltzmann method to the one-dimensional modeling of pulse waves in elastic vessels

The one-dimensional nonlinear equations for the blood flow motion in distensible vessels are considered using the kinetic approach. It is shown that the Lattice Boltzmann (LB) model for non-ideal gas is asymptotically equivalent to the blood flow equations for compliant vessels at the limit of low Knudsen numbers. The equations of state for non-ideal gas are transformed to the pressure-luminal area response. This property allows to model arbitrary pressure-luminal area relations. Several test problems are considered: the propagation of a sole nonlinear wave in an elastic vessel, the propagation of a pulse wave in a vessel with varying mechanical properties (artery stiffening) and in an artery bifurcation, in the last problem Resistor–Capacitor–Resistor (RCR) boundary conditions are considered. The comparison with the previous results shows a good precision.     

超声搅拌磁流变抛光液的声场仿真分析

为了探究超声搅拌磁流变抛光液的制备及优化工艺,利用多物理场数值计算方法,建立了超声搅拌磁流变抛光液的声场仿真模型,并进行了频域分析。研究了不同液位深度、超声变幅杆探入深度,不同功率下磁流变抛光液的声场分布。通过测量磁流变抛光液的声场强度对声场仿真进行了验证。结果表明：随着距变幅杆距离的增加,声强逐渐减弱,高声强区域主要分布在换能器轴线附近。声强在距变幅杆20mm范围内急剧衰减,变幅杆最佳探入深度为10mm,增大功率有助于空化区域的扩大。声场仿真结果与实验测量结果基本一致,对磁流变抛光液的制备提供了数值计算基础。     

Scattering coefficients for a sphere in a visco-acoustic medium for arbitrary partial wave order

Analytical solutions are reported for the scattering coefficients of a solid elastic sphere suspended in a viscous fluid for arbitrary partial wave order. Expressions are derived for incident compressional and shear wave modes, taking into account the viscosity of the surrounding fluid and resultant wave mode conversion. The long compressional wavelength limit is employed to simplify the derivation, whereas no restriction is placed on the shear wavelength in the fluid compared to the particle dimension. The analytical approximations are compared with numerical results obtained from matrix inversion of the boundary equations and agree within the validity domain of the solutions.     

A Runge-Kutta–based WENO reconstruction on moving mesh for compressible fluids

In this paper, we construct a high-order moving mesh method based on total variation diminishing Runge-Kutta and weighted essential nonoscillatory reconstruction for compressible fluid system. Beginning with the integral form of fluid system, we get the semidiscrete system with an arbitrary mesh velocity. We use weighted essential nonoscillatory reconstruction to get the space accuracy on moving meshes, and the time accuracy is obtained by modified Runge-Kutta method; the mesh velocity is determined by moving mesh method. One- and two-dimensional numerical examples are presented to demonstrate the efficient and accurate performance of the scheme.